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1.  INTRODUCTION 


Long  path  atmospheric  transmission  in  atmospheric  window  bands  such  as 
the  3-5  vun»  8-12  ym,  and  millimeter  wave  regions  can  be  limited  by  absorption 
in  the  line  wings.  There  are  now  a  number  of  critical  military  applications, 
including  infrared  search  and  track  detection  of  cruise  missiles  and  bombers, 
infrared  imaging  from  space,  millimeter  wave  communications ,  laser  target 
designation,  and  laser  weapons  beam  propagation  which  require  detailed 
knowledge  of  absorption  lineshapes  and  their  influence  on  atmospheric 
radiation  transport.  There  are  also  potential  geophysical  applications  such 
as  atmospheric  sounding  and  trace  species  monitoring  which  require  better 
lineshape  knowledge  than  that  currently  available.  Lastly,  there  is  still 
active  interest  in  energy  fluxes  in  planetary  atmospheres  and  the  global 
radiation  budget,  which  motivated  the  first  studies  of  atmospheric  window 
transmission. 

These  needs  have  resulted  in  the  development  of  a  large  body  of 
atmospheric  transmission  calculation  technology,  primarily  through  the 
Air  Force  Geophysics  Laboratory  (AFGL).  Development  continues  on 
sophisticated  computer  codes  such  as  the  AFGL  LOWTRAN  family  for  low  spectral 
resolution  and  FASCODE,  based  on  the  AFGL  line  parameters  compilation.  For 
both  codes,  the  largest  uncertainties  are  in  continuum  components  which 
dominate  in  window  bands. *  This  has  generated  continued  interest,  including 
diode  laser  and  conventional  absorption  spectroscopy  experiments,  AFGL 
workshops  on  high  resolution  atmospheric  transmission  and  yearly  AFGL  Symposia 
on  Seeing  Through  the  Atmosphere,  as  well  as  an  International  Workshop  on 
Atmospheric  Water  Vapor. ^ 


This  is  the  tenth  and  final  R&D  Status  Report  for  Contract  Number 
F19628-85-C-0183,  "Molecular  Line  Shape  Effects  on  Atmospheric  Window 
Absorption".  It  covers  work  through  contract  termination  on  February  3,  1988, 
and  includes  a  review  of  all  work  on  the  contract  since  the  start  date  of 
September  9,  1985.  The  objective  of  the  program  was  to  develop  theoretical 
methods  for  predicting  the  shape  of  molecular  absorption  lines  in  the  wing 
regions.  The  investigation  included  but  was  not  limited  to  asymptotically 
correct  forms  valid  in  the  far  wings  and  near  line  center  (to  give  line 
widths).  Its  motivation  was  to  assess  the  effect  of  line  shapes  on 
atmospheric  window  absorption,  so  that  systems  of  particular  interest  included 
water-water,  water-nitrogen,  and  carbon  dioxide-nitrogen.  Particular  emphasis 
was  placed  on  the  temperature  dependence  of  line  shapes  and  widths,  since 
these  are  among  the  most  sensitive  experimental  tests  of  detailed  theories. 


1. 1  Background  and  Research  Objectives 

The  shape  of  infrared  spectral  line  near  line  center  is  known  to  be  well 
described  by  the  Lorentz  function.  The  problem  of  calculating  absorption  near 
the  center  of  a  strong  line  is  therefore  that  of  computation  of  the  Lorentz 
width.  Microscopic  theoretical  formulations  are  available  which  allow 
evaluation  of  a  width  operator  in  terms  of  intermolecular  potentials.  From 
the  assumptions  used  in  deriving  easily  evaluated  expressions  for  the  width 
operator  which  are  valid  near  line  center,  (or  simply  from  the  requirement 
that  the  integrated  absorption  over  the  entire  lineshape  must  be  finite),  it 
is  clear  that  at  some  point  in  the  line  wing,  the  Lorentz  function  can  no 
longer  be  correct.  The  problem  of  predicting  the  wing  functional  form  as  well 
as  absolute  magnitude  has  attracted  much  attention,  but  work  to  date  has 
produced  either  formal  derivations  which  stop  short  of  numerical  results,  or 
calculations  of  wing  behavior  which  are  highly  approximate. 
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Davies,  Tipping  and  Clough^»4  (DTC)  developed  an  explicitly  time 
dependent  theory  for  spectral  line  broadening  based  on  the  dipole  moment 
autocorrelation  function.  In  contrast  to  earlier  treatments,  the  DTC  theory 
is  fully  quantum  mechanical  and  rigorously  satisfies  the  fluctuation 
dissipation  theorem  on  a  microscopic  level.  In  addition,  whereas  alternate 
quantum  mechanical  approaches^ » 6  (e.g.  those  based  on  close-coupling 
calculations  or  various  approximations  thereof)  are  typically  limited  to 
line  widths  for  smaller  molecular  systems,  the  DTC  formulation  provides  the 
entire  frequency  dependent  absorption  profile  for  more  complicated  molecular 


systems . 


The  DTC  formulation  relies  on  two  basic  approximations.  First,  it  is 
assumed  that  the  many  body  collisional  problem  for  far  wing  absorption  can  be, 
at  least  qualitatively,  reduced  to  treatment  of  the  scattering  interaction 
between  a  single  pair  of  absorbing  (radiator)  and  perturbing  molecules  and 
then  simply  multiplying  by  the  total  number  of  such  pairs.  Second,  it  is 
assumed  that  the  anisotropic  intermolecular  interaction  is  sufficiently  weak 
to  permit  the  time  dependent  scattering  transition  probabilities  to  be  treated 
using  pertubation  theory.  Although  there  is  some  evidence  that  the  so-called 
single-perturber  approximation  does  qualitatively  reproduce  observed  spectral 
features  at  low  gas  densities , ^ >8  the  reliance  on  perturbation  theory  is  more 
problematic . 


One  unresolved  issue  associated  with  perturbative  treatments  of  line 
broadening  is  the  implicit  restrictions  on  the  intermolecular  interactions 
which  may  contribute  to  far-wing  absorption.  For  example,  scattering 
resonances  and/or  complex  formation  cannot  be  adequately  treated  within  the 
DTC  formulation;  rather,  it  is  generally  assumed  that  the  dominant 
interactions  responsible  for  far-wing  absorption  arise  from  long  range 
anisotropic  interactions.  However,  even  in  this  case,  perturbation  theory 
breaks  down  as  the  impact  parameter  approaches  zero.  Practical  applications. 


therefore,  require  ad-hoc  cutoff  functions  which  ensure  physically  reasonable 
values  for  the  absorption  coefficient  in  strong  collision  regimes. 

Davies  and  Fahey  used  the  DTC  theory,  assuming  dipole-dipole  and  dipole- 
quadrupole  interactions  and  neglecting  line  coupling  effects,  to  treat 
far-wing  continuum  absorption  for  the  H2O-H2O  system.^  We  have  also  used  DTC 
theory  to  study  this  system.  In  both  cases  it  was  found  that  the  calculated 
profile  failed  to  agree  with  available  experimental  data.  In  particular,  the 
theoretical  far-wing  absorption  cross  sections  had  the  wrong  temperature 
dependence.  In  addition,  our  analysis  indicated  that  the  cutoff  functions 
needed  to  ensure  well  behaved  absorption  profiles  under  strong  intermolecular 
interactions  were  affecting  the  H2O-H2O  line  shapes  out  to  1000  cm-*  off  line 
center. 

On  the  basis  of  this  analysis  it  was  suggested  that  additional  research 
was  needed  to  develop  theoretical  methods  which  were  either  nonperturbative 
or,  at  least,  less  dependent  upon  ad  hoc  cutoff  functions  to  deal  with  strong 
collisions  between  the  radiating  and  perturbing  species.  In  addition,  it  was 
also  suggested  that  the  theoretical  methods  and  corresponding  computer  codes 
should  be  more  general  with  respect  to  the  type  of  intermolecular  interactions 
with  which  they  can  deal. 

Accordingly,  our  research  has  focused  on  developing  two  theoretical 
methods  which  showed  promise  for  achieving  these  goals.  Both  methods  are 
based  on  the  dipole  autocorrelation  function  formulation  and  assume  uncoupled 
lines  and  the  binary  collision  approximation.  The  first  method,  discussed  in 
Section  4,  is  based  on  a  generalization  of  the  leavitt-Korff  cutoff-free 
impact  theory  of  line  broadening.^  This  method,  although  based  on 
perturbation  theory,  uses  a  linked-cluster  resummation  of  the  perturbation 
expansion  to  yield  absorption  cross  sections  which  remain  well-behaved  in 
strong  collision  regimes  without  resorting  to  ad  hoc  cutoff  functions.  The 


second  technique,  the  Recursive-Residue-Generation  Method  (RRGM)IO’H  is  a 
nonperturbative  method  which  can  treat  arbitrary  intermolecular  interactions. 

1.2  Report  Outline 

Section  2  discusses  the  basic  approach  and  common  approximations  invoked 
by  the  theoretical  methods  described  in  this  report.  Section  3  describes  the 
application  of  the  DTC  method  (developed  in  Appendix  A)  to  H£0  self-broadening. 
Section  4  describes  the  GLK  method  and  gives  results  for  the  H2O 
self-broadening.  Section  5  describes  the  RRGM  method  which  is  illustrated  by 
an  application  to  HCl-Ar.  A  summary  is  presented  in  Section  6. 

Some  mathematical  details  are  treated  in  the  Appendices.  In  addition  to 
a  formal  description  of  the  DTC  lineshape  in  Appendix  A,  Appendix  B  and  C 
contain  a  description  of  the  matrix  elements  of  high  order  multipole 
interaction  and  computational  techniques  for  radial  matrix  elements, 
respectively. 
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2.0  OVERVIEW  OF  THEORETICAL  METHODS 


2.1  Autocorrelation  Function  Formulation  of  t 


lLrgfllM 


The  absorption  coefficient  per  unit  path  length,  a(w),  for  a  dilute  gas 
is  related  to  the  thermally  averaged  dipole  autocorrelation  function  for  the 
system,  $(t),  by1^ 


a(u)  = 


(l-e"eRw) 


e"iut  $(t) 
2tiR  6 


(2-1) 


where  3  =  1/kgT,  nr  is  the  number  density  of  radiating  dipoles,  and  $(t)  is 
the  dipole  autocorrelation  function  for  a  single  radiator  interacting  with  Np 
perturbers , 


$(t)  =  Tr[p(0)y  .  fe 


iHt/fi 


-iHt/R, 


(2-2) 


where  p(0)  is  the  initial  density  matrix,  y  is  the  initial  dipole  moment 
operator  and  H  is  the  Hamiltonian  operator  for  the  system.  We  shall  assume  a 
canonical  system  so  that 


P ( 0 )  =  p(H)  = 


Tree’6") 


(2-3) 


is  the  canonical  density  operator  with  partition  function  Qp.  In  addition, 
the  present  report  will  specifically  treat  foreign  broadening  by  perturbing 
molecules  which  do  not  absorb  in  the  spectral  regions  of  interest.  The  dipole 
moment  operator  is  then  y  =  yr  where  yr  is  the  dipole  moment  of  the  absorbing 
species  (radiator).  However,  self-broadening  may  be  treated  in  basically  the 
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same  fashion. ^  Lastly,  defining  the  complex  time  variable  z  =  t-i0R  and  using 
the  cyclic  invariance  of  the  trace  operation,  Eq.  (2.2)  can  be  rewritten  as 


T  (  iHt/R  ■+  -iHz/R,  . 

$(t)  =  Tr [p  .  {  e  p  e  }  ] 

P 


(2-4) 


where  the  Boltzmann  statistical  weight  has  been  incorporated  into  the  time 
evolution  of  the  dipole  moment  operator. 


All  calculations  in  this  report  will  assume  uncoupled  lines  and  the 
binary  collision  approximation.  The  uncoupled  line  approximation  assumes  that 
there  are  not  collision  induced  correlations  between  zero-order  transitions. 
Then,  the  autocorrelation  function  may  be  expressed  as  a  sum  of  correlation 
functions  for  each  dipole  transition,  i.e., 

Ht)=l  l  i  At) 

i  f  1 

where  i  and  f  refer  to  initial  and  final  states.  The  binary  collision 
approximation  assumes  that  the  Np  perturber  autocorrelation  function  is  simply 
the  product  of  Np  identical  single-perturber  autocorrelation  functions.  In 
this  case  it  is  shown  in  Appendix  A  that 

iw.fT  -nflT  (t) 

$if (t)  a  e  e 

where  w^f  is  the  transition  frequency,  np  the  perturber  intensity,  fl  the 
system  volume,  and  is  approximately  a  properly  normalized  single  perturber 
correlation  function. 
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3.0  APPLICATION  OF  DTC  THEORY  TO  H20  SELF-BROADENING 


3.1  H20  Absorption  In  the  10  pm  Region 

For  atmospheric  "window"  transmission,  H20  absorption  is  particularly 
important  in  the  10  pm  region.  It  is  generally  thought  that  absorption  in 
this  region  is  due  to  wings  of  lines  in  the  6.3  pm  and  pure  rotation  bands. 

It  has  also  been  speculated  that  dimers  (H20  •  H20  or  H20  •  N2)  or  larger 
clusters  may  be  responsible  for  much  of  the  absorption.  Since  line  wing  and 
dimer  absorption  have  the  same  pressure  dependence,  and  indeed,  at  observed 
temperatures  could  have  the  same  continuum  spectral  appearance,  it  is  not 
possible  to  distinguish  between  different  mechanisms  from  present  experimental 
or  theoretical  work. 

Line  wing  absorption  is  often  reported  as  the  ratio  of  observed 
absorption  to  that  expected  for  a  Lorentzian  lineshape.  This  ratio  is  denoted 
X.  At  this  time,  the  general  observation  is  that  water  continuum  absorption 
in  the  extreme  wings  lies  below  the  absorption  predicted  by  summing  the 
Lorentzian  line  wings  from  pure  rotational  transitions.^^  This  is  illustrated 
in  Figure  3.1  which  displays  the  temperature  dependent  X  factor  corresponding 
to  the  ratio  of  the  experimental  absorption  to  the  absorption  computed  by 
summing  Lorentzians  for  all  distant  lines  on  the  HITRAN  tape,  after 
subtracting  out  all  nearby  lines. 1^.15 

In  Figure  3.1,  the  experimental  X  factors  at  1000  cm”^  range  from  0.1  to 
0.6  and  decrease  with  increasing  temperature.  To  understand  the  theoretically 
predicted  temperature  dependence  of  far  wing  absorption,  the  DTC  formulation 
was  used  to  treat  a  few  isolated  H20  rotational  lines.  The  model  system  and 
results  are  presented  in  Section  3.2.  The  predicted  temperature  dependence  is 
discussed  in  Section  3.3. 
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Figure  3.1.  Experimental  Water  Continuum  As  Lorentzian  Correction, 
at  1000  cm"* 
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The  DTC  method,  assuming  uncoupled  lines  and  the  binary  collision 
approximation,  was  used  to  treat  a  few  isolated  lines  for  the  H2O/H2O  system. 
The  zeroth-order  Hamiltonian  is  given  by 


Ho  =  <  +  +  «RP 

where  H°,  H°  are  rigid  rotor  Hamiltonians  for  the  H2O  molecules  with 
eigenstates  |  j,T,m>  and  eigenvalues  £(j,T).  The  operator  H°R  describing  the 
relative  motion  of  the  H2O  molecules  is  given  by 

Ho  =  l_  R2  a_  lL_ 

HRp-2yR23RR  3R  2yR2 

where  y  is  the  reduced  mass  for  two  point  particles  with  reduced  mass  y  =  H2O 
R  is  the  intermolecular  separation,  and  L  is  the  total  angular  momentum. 

Using  the  spherical  harmonics  satisfying 


L2!  £m>  =  R2£(£  +  1)  |  £m> 
L  |  £m>  =  Rm  |  £m> 


the  stationary  states  of  H°p  can  be  expressed  as  |  fcy£>  = 
|  f">  i  £$>  where  the  radial  eigenfunctions  satisfy 


l  -  li  +  +  VR>1  I  «>  -  <1  «> 

2yR  2yR 


For  the  present  calculation  the  isotopic  intermolecular  interaction  has  been 
modeled  by  a  6-12  Lennard- Jones  potential 


Vo  (R)  -  Vo  [  (6/R) 12  -2(6/R)6] 

The  anisotopic  coupling  interaction  responsible  for  modifying  the  H2O 
rotational  lines  was  represented  as  an  multiple  expansion  which  includes  terms 
through  electrostatic  dipole/quadrupole. 

v(R,e,*)  =  vDD(R,e,*)  +  vDQ(R,6,0)  +  vqdcr,g,*) 

In  a  cartesian  coordinate  system  these  terms  are  given  by 


VDD  =  h  ^1  *  *2  •  3(h  ‘  R)(^2  *  R)I 
K 

VDQ  =  h  [2r»Q2«'u1  -  5  (R.y1)(R  Q^R)] 
2R 


where  yj,  yj  and  Qj  and  Q2  are  dipole  and  quadrupole  cartesian  tensors  for 
the  radiator  (1)  and  perturber  (2),  respectively.  Parameter  values  are  given 
in  Reference  A. 

The  DTC  chi-factors  at  1000  cm"l  for  the  (A,  1,5)  -+  (5,2,5)  H2O  rotational 
transition  are  shown  in  Figure  3.2.  In  addition  to  being  an  order  of 
magnitude  smaller  than  the  experimental  results  the  theoretical  chi-factors 
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also  show  a  positive  temperature  dependence.  We  were  able  to  show  that  this 
is  not  unexpected,  using  simple  qualitative  arguments.  We  review  that  work 
here. 
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Simple  models  to  correlate  DTC  results  can  be  found  in  analytic 
approximations  to  Anderson  theory  for  line  widths.  A  particularly  useful 
review  is  given  by  Birnbaum  in  Volume  12  of  Advances  in  Chemical  Physics.  It 
is  shown  that  the  line  broadening  cross  section  consists  of  two  factors.  The 
first  is  a  coefficient  which  is  proportional  to  a  power  of  the  temperature. 

The  second  is  a  function  of  a  resonance  factor  which  is  proportional  to  the 
total  change  in  rotational  energy  of  the  absorbing  and  colliding  molecules. 

Much  rule-of  thumb  analysis  of  linewidth  temperature  behavior  is  done  by 
assuming  that  in  the  important  line  broadening  collisions,  total  changes  in 
rotational  energy  and  therefore  the  resonance  factor  are  small,  and  its 
function  remains  near  one.  In  that  case,  the  power  of  the  temperature 
corresponding  to  the  DTC  definition  of  the  linewidth  is  (n-3)/(n-l),  where  n 
is  the  exponent  varies  from  zero  to  one-half  as  the  potential  is  varied  from 
dipole-dipole  (n=3)  to  hard-sphere.  The  DTC  predictions  on  resonance  lie 
within  this  range,  but  for  non-resonant  photons  the  predictions  increase  much 
more  rapidly  with  increasing  temperature. 

In  fact  there  is  no  reason  for  the  water-water  system  to  follow  the  above 
simple  dependence.  For  non-resonant  photons,  the  resonance  factor  (which  is 
essentially  the  collisional  change  in  momentum  required  to  conserve  energy) 
contains  the  difference  between  the  photon  and  transition  energies.  Large 
contributions  are  obtained  with  transitions  with  small  resonance  factors, 
except  that  each  transition  is  weighted  by  the  Boltzmann  factor  for  the 
initial  rotational  state  of  the  perturber.  That  weighting  determines  an 
average,  non-zero  resonance  factor. 
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By  taking  the  resonance  factor  function  to  be  an  exponential  of  a  power 
of  the  inverse  temperature,  one  can  deduce  averages  of  the  translational 
energy  change  during  a  collision  by  fitting  the  resulting  exponential 
expression  to  the  predicted  DTC  temperature  dependence  for  the  (4,1,4)  - 
(5,0,5)  transition.  The  variation  in  this  quantity  with  varying  photon 
energy,  with  excluding  or  including  dipole-quadrupole  forces,  and  between  the 
initial-  and  final-state  contributions  is  with  only  one  exception  in 
accordance  with  the  expectation  that  larger  average  translational  energy 
changes  will  be  associated  with  smaller,  more-temperature  dependent  terms. 

The  one  possible  discrepancy,  between  the  two  potential  choices  for  the 
initial  state  term  and  a  1000  cm"*  photon,  may  be  due  to  the  different 
weighting  of  perturber  rotational  transitions  in  the  two  cases. 

In  order  to  investigate  the  details  of  which  collisional  transitions  are 
most  important,  we  developed  a  computer  program  (S2TERM)  which  read  the  same 
input  files  as  the  DTC  HWIDTH  program,  for  a  single  translational  energy  and 
angular  momentum  quantum  number  (£).  S2TERM  examines  the  contributions  to  the 
linewidth  function  from  all  rotational  transitions  in  absorber  and  pertuber, 
then  prints  out  data  for  only  those  whose  contribution  is  above  a  certain 
fraction  of  the  total.  The  simplification  is  impressive  -  for  example,  for 
the  (11,3,8)  -  (12,4,9)  dipole-quadrupole  and  quadrupole-dipole  calculations 
for  H  =  0,  400  cm”l  translational  energy  and  a  1211.3  cm-*  photon  energy,  47 
and  59  terms  out  of  10,000  give  over  half  and  almost  two-thirds  of  the  totals. 

The  data  printed  out  in  the  S2TERM  lists  include  the  changes  in  absorber 
and  perturber  rotational  energy  and  relative  translational  energy,  the 
Boltzmann  factor  for  the  perturber  initial  rotational  state,  the  transition 
strengths  for  both  rotational  changes,  and  the  matrix  element  appropriate  to 
the  residual  change  in  translational  energy.  It  is  the  product  of  the  last 
four  terms  that  determines  the  contribution  from  a  given  pair  of  rotational 
transitions . 


The  sum  of  the  three  collisional  energy  changes  satisfy  an  energy 
conservation  relation  with  the  difference  between  the  photon  and  transition 
energies.  The  relation  is  different  for  the  initial  and  final  state  terms. 
Using  the  notation  in  the  GTE  final  report  is 

(tj  -  Ej,  )  +  (  Ej  -  Ejt)  +  (Ek  -  Ek.)  -  (w  -  wfi)  =  0 

for  the  initial  state  term,  and 

(E j  -  Ej  )  +  (Ej  -  Ej.)  +  (ER  -  Ek.)  +  («  -  Wfi)  *  0 

for  the  final  state.  If  we  designate  the  first  three  terms  as  decon,  these 
relationships  can  be  illustrated  as  follows. 
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It  can  be  seen  that  for  initial  absorber  states  of  relatively  low  energy 
much  of  the  downward  change  in  system  energy  required  to  accommodate  the 
off-resonance  photon  will  have  to  come  either  from  perturber  transitions  from 
high  rotational  states  with  small  Boltzmann  weightings,  or  from  large  changes 
in  translational  energy  giving  rise  to  small  matrix  elements.  There  are  no 
such  constraints  on  the  final  state  term,  so  that  in  the  far  (blue)  wing  the 
important  terms  in  its  sum  are  one  to  two  orders  of  magnitude  larger  than 
those  contributing  to  the  initial  state  cross  section.  However,  the  energy 
needed  to  carry  the  state  f  into  f'  should  also  appear  in  a  Boltzmann 
weighting,  and  the  inverse  of  the  factor  e[-BR(w-Uf^) ]  (a  detailed  balance 
factor)  multiplying  the  initial  state  term  restores  the  proper  ratio. 


In  the  low-J  (strong  line)  cases  studied,  the  final  and  detailed  balance 
weighted  initial  state  terms  are  of  the  same  order  of  magnitude.  Although  the 
detailed  balance  factor  has  a  very  strong  negative  temperature  dependence, 
much  like  that  actually  observed  in  water-water  far  wing  absorption,  the 
initial  state  term  for  low  J's  has  a  very  strong  positive  temperature 
dependence  as  well  as  a  small  magnitude  (both  due  to  perturber  rotational 
Boltzmann  factors).  Meanwhile,  the  contributions  from  the  final  state  term 
come  mostly  from  higher  J  transitions  as  well,  since  high  transition  strengths 
and  large  translational  matrix  elements  (when  a  transition  between  two  high 
J's  can  mal  up  most  of  the  energy  difference)  make  up  for  smaller  perturber 
rotational  populations.  The  result  of  these  trends  is  the  strong  positive 
temperature  dependence  seen  in  the  GTE  predictions. 

This  temperature  dependence  is  less  strong  when  dipole-quadrupole  forces 
are  included.  It  is  reasonable  to  expect  some  further  improvement  is 
quadrupole-quadrupole  forces  were  included.  It  is  also  less  strong  for 
absorbers  in  high  initial  rotational  states.  S2TERM  analysis  shows  these 
trends  to  be  explained  by  the  above  discussion,  specifically  by  changes  in  the 
perturber  J  states  which  can  contribute  and  the  resulting  changes  in  cross 
section  magnitude  and  temperature  dependence.  We  ran  GTE  far  wing  predictions 
for  two  high-J  and  two  low-J  transitions,  and  found  no  behavior  which  could 
not  be  correlated  as  described  above.  Although  high-J  lines  have  negative 
temperature  dependence  and  super-Lorentzian  wings  as  do  the  observed  wings, 
their  intensities  are  much  too  small  for  them  to  make  significant 
contributions . 

One  conclusion  which  can  be  drawn  from  the  above  investigations  and  our 
earlier  work  on  the  unimportance  of  bound  and  quasi-bound  states  is  that  the 
only  source  of  strong  negative  temperature  dependence  of  wing  absorption  is 
the  detailed  balance  factor.  If  this  were  true,  one  would  expect  that  the 
observed  temperature  dependence  would  become  more  strongly  negative  with 
increasing  frequency  difference  from  line  center.  Our  analysis  of  the  data 


reported  by  Burch  seems  indeed  to  show  a  variation  in  the  temperature  exponent 
that  parallels  that  obtained  from  fitting  the  exponential  in  inverse 
temperature  to  a  T"n  form. 

The  observations  are  not  so  precise  that  any  agreement  would  constitute 
proof,  and  in  fact  the  observed  exponent  seems  to  be  about  one  higher  than 
that  obtained  from  the  detailed  balance  term  alone,  without  any  modifications 
by  linewidth  factors  which  still  will  increase  with  temperature.  On  the  other 
hand,  comparison  with  experiment  does  not  seem  to  rule  out  the  possibility 
that  a  predictive  method  based  on  the  GTE  model,  but  with  additional 
transitions  included  to  improve  linewidth  term  magnitudes  and  tempei  .ture 
dependences,  might  in  fact  be  an  accurate  description  of  wing  absorption. 


4.0  CUT-OFF  FREE  THEORY  FOR  ROTATIONAL  LINESHAPES 
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4. 1  Generalized  Leavitt-Korf f  Theory  of  Line  Broadenii 


Leavitt  and  Korff^  derived  a  variation  on  Anderson's  impact  theory1^  of 
line  broadening  which  remains  well-behaved  as  the  impact  parameter  approaches 
zero.  Consequently,  applications  do  not  require  ad  hoc  cutoff  functions  to 
preserve  unitarity  in  strong  collision  regimes.  There  are  two  key  features  to 
their  derivation.  First,  it  is  assumed  that  the  time  evolution  operator  is 
both  energy  conserving  and  invariant  with  respect  to  orientations  of  the 
radiating  and  perturbing  molecules.  Second,  a  linked-cluster  theorem 
for  degenerate  states  is  used  to  decompose  each  order  in  the  perturbation 
series  in  terms  of  its  linked-cluster  expansion.  Subsequent  factorization  and 
resummation  then  yields  an  expression  for  the  interruption  function  in  which 
Anderson's  S2~functions ^  have  been  exponentiated.  In  this  section,  the 
Leavitt-Korff  half-width  theory  is  generalized  for  frequency  dependent 
lineshapes . 

In  Equation  (A. 14),  the  Hamiltonian  H0  (for  convenience  we  will  drop  the 
superscript  and  subscript  s  notation)  for  a  single  perturber  interacting 
with  the  radiating  molecule  is 


Ho  =  H?  +  I& 


(4-1) 


where  H  is  the  rigid  rotor  Hamiltonian  for  the  perturber  with  eigenstates 
P 

H“  I  XF  >  =  H°  |  jp  Tp  mp  >  =  E.^  |  jp  ,p  mp  >  (4-2) 

is  the  Hamiltonian  operator  for  two  structureless  particles  interacting 
through  a  two-body  local  isotropic  potential  V°(R) 


- 1-5;  ?  2  +  v°«,] 


The  stationary  states  of  H  are 

RP 


(4-3) 


I  S  >  =  “ip  1 1  ■(  ‘  > 


(4-4) 


Using  these  zeroth-order  states,  the  trace  in  Equation  (A- 14)  can  be  evaluated 
to  give 

Tr  [ . . .  ]  =5;  III  p(E.  t  +  e)  <  j.t.m'  \Z  |  U(z)  |  j^Vp  ^ 


Xp  k  X’  k’ 


x  <jfTf®fxp  I  u+(t)  I  xp^  *  (4“5) 


As  a  first  approximation  to  Equation  (4.6a)  it  is  assumed  that  the  matrix 
elements  of  the  time  evolution  operators  are  approximately  given  by 


<jiTimi  XpS  I  0(*>  ^  6XpX-  6^,  5i  <z> 


<^fTfmf  XP  I  U+(z)  ^fTfmfXP  ^  >  ~ 


«  5X  x.  6  .  0  (t) 

Vf  XPXP  k,k’  f 


(4-6a) 


(4-6b) 


where 


^(z)  »  f(2j1+l)(2jp+l)(2«  +  l)]'1  l  l  l  <  jjT^Xp  t  |  U(z)  |  jjt^Xp  it  > 

“i  nP  mt  (4-7a) 


(Z)  =  [(2jf+l)(2jp+l)(2£  +  l)]"1  l  l  l  <  j^fXp  i  I  utt)  |  £  > 

”f  mP  (4-7b) 

As  in  Reference  (1),  the  basis  of  this  approximation  is  that  the  resulting 
expression  for  autocorrelation  function  reduces  to  that  given  by  DTC  when  the 
perturbation  is  weak.  Substituting  Equation  (4-6)  into  Equation  (4-5)  gives 

Tr[...J  =  6  ,  6  t  ll  l  f  de  (2j+l)(2£+l)  p(E  +  e)  U  (z)  0+(t) 

mimi  “f^f  jp  tp  l  J  P  JPTP  1  f 

Then,  substituting  Equation  (4-8)  back  into  Equation  (A-14)  and  performaing 
the  summations  of  magnetic  quantum  numbers  gives 


Q.f(t)  =  l  l  l  f  dt  (2j  +l)(2t+l)  p(E  +  e)  U  (z)  U+(t)  (4-9) 

11  ;  ^  o  r  Jplp  A  * 


iP  t 


The  matrix  elements  and  Uf  can  be  evaluated  using  the  linked-c luster  theorem 
for  degenerate  states  as  described  in  Reference  (9)  giving 


V  S  .(z) 
L  n,i 


Ui(z)  =  e 


(4-10a) 


l  Sn  #<*> 

u  n,  f 


(t)  =  e" 


where  the  Sn  functions  are 


n  r 

Sn,i(z)  =  C-i/K)  /  dZl  . J  dzn  [(2j.+  l)(2Jp+  1)(2£  +  1)] 


-1 


(4-10b) 


(4-lla) 


xl  I  I  <  X.XpS  |  [V(Z1)  ...  V(zn))^0)  |  X.Xpk 


m  m  mQ 
i  p  t 


Sn,f(t)  =  dtl  •••/  dtn  [(2jf+  1)(2jp+  1)(2£  +  1)1_1 

o  o 

X  l  l  l  <  XfXplc  !  [V(ti)  ...  V(tn)]£0)  I  XfXpk  >  (4-llb) 

mf  “p  ”£ 

As  in  Reference  (9),  the  subscript  L  indicates  that  in  evaluating  the  quantity 
in  brackets  all  subproducts  which  are  diagonal  in  all  quantum  numbers  and 
rotationally  invariant  in  the  spaces  of  the  radiator,  perturber  and  radial 
vector  1$  are  excluded,  while  the  superscript  (0)  indicates  that  the  entire 
quantity  is  rotationally  invariant.  As  noted  in  Reference  (9),  the  first  order 
terms  Sj^  and  S^f  are  complex  and  only  contribute  to  the  line  shift.  The 
S£~functions  are  then  the  first  terms  to  contribute  to  line  broadening.  These 
are  given  by 


S2,i<Z>  = 


-(l/K)2/dZj  f  dz2[(2ji+  l)(2jp+  1)(2£+1)]’1  x  ^  ^ 
o  o  Xi  XP  k' 

x  l  l  l  <  X.Xp£  |  V(Zl)  |  Xj  X*  £’><Xj  X*  k'  |  V(z2)  |  X^S  > 

«i  »p  mt 


(4-12a) 


S2,f(t>  “ 


-d/B) 


f  dti  f1  dtz[(2jf+  1) (2 jp4-  1)(2£+1)]'1  x  l  l  l 

o  o  xi  x:  P 


mfjmm 


X  l  l  l  <  X.Xpk  I  V(tx)  |  x:  XJ  k'xx:  X*  k-  I  V(t2)  I  X.Xpk  > 


mi  mp  m£ 


(4-12b) 


where  the  primes  on  the  summations  indicate  that  no  diagonal  terms  are  to  be 
included  in  the  summation.  Then,  up  through  second  order  in  the  interaction, 
the  single-perturber  autocorrelation  function  is 


=  I  l  l  J  de  P<Ei  t 


+  8)  e 


S2,i(z)  +  S2,f(t> 


(4-13) 


Jp  tp  £ 


Details  for  evaluating  matrix  elements  of  the  anisotropic  interaction  potential 
V,  assuming  a  multipole  expansion  are  given  in  Appendices  B  and  C,  for 
arbitrary  orders  of  the  multipolar  tensors. 


The  generalized-Leaviff-Korff  approach  was  used  to  treat  H2O 
self-broadening  using  the  molecular  Hamiltonian  given  in  Section  3.2. 

Figure  4.1  shows  chi-factors  obtained  from  the  DTC  and  GLK  methods  out  to  900 
cm“l  off  line  center,  assuming  dipole-dipole  interactions  only.  At  100  cm”^ 
the  DTC  and  GLK  results  agree.  For  higher  frequencies,  however,  the  GLK 
chi-factors  are  2-3  times  smaller  than  those  given  by  DTC. 


5.0  RECURSIVE  RESIDUE  GENERATION  METHOD 
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5 . 1  General  Description  of  the  Method 

RRGM  is  a  nonperturbative  method  for  constructing  time  dependent 
transition  amplitudes  that  is  particularly  well-suited  to  problems  with  large, 
but  sparse  matrices . 1° « H  In  this  method,  the  Hamiltonian  for  the  system  is 
first  tridiagonalized  using  the  Lanczos  Recursion  Algorithm. ^  Subsequent 
diagonalization  of  the  tridiagonal  (Jacobi)  matrix  yields  the  residues  and 
eigenvalues  needed  to  determine  a  given  transition  amplitude. 

If  {  |  a>|  and  {Ea}  are  the  eigenstates  and  eigenvalues  of  the 
Hamiltonian  H,  the  time  dependent  transition  amplitude  between  two 
nonstationary  states  |  i>  and  |  f>  is 

Afi(t)  =  <f  I  e”lHt/R  |  i>  =  l  <f  |  «>«x  |  i>e"lEat/R=  £  Rfi(a)  e"iEat/R  (5-1) 

a  a 

where  Rfi(«)  =  <f  |  a  >  <  a  |  i>  is  the  off  diagonal  residue  corresponding  to  a 
simple  pole  Ea  for  the  Green's  function, 


;fi(a)  -  <£  I  (z-H)”1  |  i>  =  l  =  l 


Rfi(a) 


(5-2) 


For  smaller  systems  the  residues  and  eigenvalues  needed  to  evaluate  Afj_(t)  can 
be  obtained  by  standard  diagonalization  methods.  For  larger  systems,  however, 
these  procedures  are  usually  prohibited  by  computer  storage  requirements. 

RRGM  circumvents  this  difficulty  because  the  actual  eigenvectors  are  never 
explicitly  generated. 


vqv; 
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It  is  easily  shown  that  upon  defining  the  vectors 


|  U()>  -  [  |  i>  +  |  f>]/V2 


I  vQ>  =  [  |  i>  -  |  f>]//2 


(5-3a 


(5-3b 


the  off-diagonal  residue  Rf^(a)  can  be  expressed  in  terms  of  the  diagonal 
residues  associated  with  the  vectors  |  uq>  and  |  vq>  as 

Rfi(ot)  =  <£  I  a><a  |  i>  =  [<uQ  |  a><a  |  u0>  -  <vQ  |  a><a  |  vQ>]/2 

=  tRu  u  <«*>  "  Rv  v  («)]/*  *  (5-4 

O’  0  O’  0 

In  addition,  if  |  ug>  is  the  first  vector  of  a  set  of  vectors  {|  uj>}  and  |  vq> 
is  the  first  vector  in  another  orthonormal  set  of  vectors  {|  vj>},  each  of 
which  tridiagonalizes  the  Hamiltonian  H,  i.e. 


U+HU  =  J 


V+HV  =  J 


(5-5a 


(5-5b 


then  the  residues  Rug,uo  and  Rvq,vq  can  be  determined  by  diagonalizing  the 
Jacobi  matrices  Jv  and  Jv.  Specifically,  if  Tu  and  Tv  are  the  matrices  which 
diagonalize  Ju  and  Jv,  respectively,  then  the  needed  residues  are  given  by 


R  (o)  »  (T  ),  (T  )  , 

u-.u.  u ^ 1 , a  u  a, 1 


R  (a)  =  (T  ),  (T  )  . 

V  V  V  1  Ct  v  u  1 


(5-6a 


(5-6b 


It  is  important  to  note  that  only  the  first  column  of  the  Tu  and  Tv  matrices 
are  needed  to  determine  all  the  residues  and  eigenvalues  required  to  compute 
the  transition  amplitude  Afj. 


5? 


Up  to  this  point  there  has  been  no  reduction  in  the  computational 
requirements.  Rather,  the  initial  problem  of  diagonalizing  the  Hamiltonian  H 
has  simply  been  exchanged  for  the  equally  formidable  problem  (for  larger 
systems)  of  determining  the  two  Jacobi  matrices.  The  key  feature  in  RRGM  is 
that  the  columns  of  the  transformation  matrices  U  and  V  (as  well  as  the 
diagonal  and  off-diagonal  elements  of  Ju  and  Jv)  are  determined  recursively 
using  the  Lanczos  tridiagonalization  algorithm^  so  that  only  two  columns  of 
the  complete  matrices  are  required  in  fast  storage.  For  example,  starting 
with  the  initial  recursion  vector  |  ug>  and  its  self-energy  ag=<ug  |  H  |  uq>,  the 
second  vector  in  the  set  is  defined  by  |  uj>  =  [(H-ag)  |  ug>]/bj  where  bj  is  the 
norm  of  the  so-called  residual  vector,  (H-ag)  [  uq>.  The  self-energy  ag  and 
norm  bj  are  the  first  diagonal  and  off-diagonal  elements  of  Ju.  Subsequent 
recursion  vectors  |  un>  and  diagonal  an  and  off-diagonal  bn  elements  of  Ju  are 
then  determined  by  the  three-term  recursion  relation 


I  u  , , >  =  [(H-a  )  j  u  >  -  bl  u  >]/b 
1  n+1  1  n  1  n  ri  n-1  J/  r 


(5-7) 


where  an  =  <un  |  H  |  un>  and  bn+i  normalizes  the  residual  vector  [(H-an)  |  un>  - 
bn  |  un_i>].  By  construction,  each  recursion  vector  is  implicitly  orthogonal 
to  all  previous  vectors.  A  second  application  of  the  Lanczos  algorithm 
starting  with  |  vQ>  as  the  initial  recursion  vector  yields  the  diagonal  and 
off-diagonal  matrix  elements  of  Jv.  Diagonalization  of  the  two  Jacobi  matrics 
is  straightforward  and  gives  the  residues  in  Equation  (5-6). 

To  p  vively  realize  the  full  potential  of  RRGM,  the  Hamiltonian  matrix 
must  be  r  tively  sparse  and/or  the  intermolecular  coupling  matrix  elements 
easily  computed  so  that  large  amounts  of  storage  are  not  required  for  coupling 
matrix  elements.  This  indeed  tends  to  be  the  case  for  the  spectral  line 
broadening  problem  for  three  reasons.  First,  the  intermolecular  interaction 
can  often  be  represented  by  an  expansion  (e.g.  multipole  expansion  or  a 
numerical  fit  to  a  Legendre  polynominal  expansion)  whose  individual  terms  can 
be  factored  into  components  which  act  only  on  the  subspaces  for  the  separated 


molecular/atomic  species  and  their  relative  motion.  Second,  selection  rules 
and/or  conservation  of  total  angular  momentum  further  decrease  the  number  of 
nonzero  coupling  matrix  elements  within  a  given  subspace.  Third,  a  portion  of 
the  interaction  coupling  resides  entirely  in  the  coupling  of  angular  momenta 
and  is  easily  computed  in  terms  of  Clebsch-Gordan  coefficients  and  6-j  symbols. 
The  net  consequence  is  that  although  the  coupling  matrix  elements  are 
significantly  more  complicated  for  line  broadening  than  in  earlier 
applications  of  RRGM,  the  basic  requirements  for  implementation  of  RRGM  should 
be  satisfied. 


5.2  Application  to  HCl-Ar 

As  a  preliminary  test  of  the  method,  RRGM  has  been  used  to  calculate  the 
dipole  autocorrelation  function  for  the  ji=5->jf=6  pure  rotational  line 
of  a  dilute  HCl-Ar  mixture.  The  HCl-Ar  system  has  been  previously  treated  by 
Boulet  and  Robert  within  a  unified  semiclassical  theory  using  curved  classical 
trajectories  and  including  all  orders  in  the  inter-molecular  interaction. 

The  results  from  this  latter  method  were  used  in  comparisons  against  the  RRGM 
results . 

5.2.1  System  Hamiltonian  and  Matrix  Elements 
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In  treating  HCl-Ar  rotational  line  broadening  we  have  used  the 
space-fixed  Arthurs  and  Dalgarno^  rigid  rotor  theory  for  scattering  of  a 
structureless  atom  by  a  X-state  rigid  rotor  diatomic.  The  coordinate  system 
in  this  approach  has  its  origin  at  the  system  center-of-mass  and  its  axes 
aligned  with  space-fixed  directions.  The  Hamiltonian,  after  separating  out 
the  center-of-mass  motion,  for  a  single  radiator  (HC1) — perturber  (Ar)  pair  is 


_2  2 

H  =  H  .  -  f-  $  +  V 

rad  2y 


(5-8) 


5-4 


where  V  is  the  Born-Oppenheimer  intermolecular  potential,  7  is  the  Laplacian 
operator  of  R  (vector  between  the  atom  and  diatomic  center-of-mass) ,  p  is  the 
reduced  mass  of  Ar  relative  to  HC1,  and  Hra(j  is  the  Hamiltonian  for  the 
nuclear  motion  of  the  radiator.  Since  we  are  interested  in  rotational  line 
broadening,  we  will  simplify  the  problem  by  neglecting  the  diatomic 
vibrational  degree  of  freedom  and  treating  HC1  as  a  rigid  rotor.  The 
eigenfunctions  of  Hra(i  then  satisfy 


H  .  Y.  =  E.  Y. 
rad  jnu  j  jm^ 


(5-9) 


where  Yjm^  is  a  spherical  harmonic  in  the  angular  degrees-of-freedom  of  r, 
(0r,0r).  Similarly,  if  L  is  the  angular  momentum  of  Ar  relative  to  HC1,  then 
the  eigenfunctions  of  the  operator  Lz  and  Lz  are  Y£,m£  (9r>$r)- 

Since  the  total  angular  momentum  J  =  Jra(j  +  L  is  conserved,  it  is 

convenient  to  couple  Y;m.  and  Y£  m£  to  form  eigenfunctions  of  Jz,  Jz, 

J  J  ’ 

JradZ  and  Lz  using  the  Clebsch-Gordan  expansion 


|  jfcJM>  =  X  l  (jm  l  m  |  jfc.M)  lim*5,  I  *■  > 

m =  -j  m  =  -£  J  J  J 


(5-10) 


In  addition,  since  RRGM  relies  of  matrix  operations,  the  translation  continuum 
will  be  discretized  by  enclosing  the  atom-diatomic  system  in  a  sphere  of 
radius  R0.  A  spherical  Bessel  basis  set, 
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is  then  used  to  span  the  translational  quasi-continuum. 


(5-11) 


In  terms  of  these  zero-order  basis  functions  the  Hamiltonian  matrix 


elements  are 


2  R2 

<j ' £ ’  j'M'n'l  H  lj*JMn>  =  6.,  [ { (Ej  +  1 “r)  fl  .  +  (5-12 
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If  the  intermolecular  potential  is  expanded  in  Legendre  polynomials 


V(R,0)  *  l  VX(R)  Px(cos9)  (5-13 


then  the  coupling  matrix  elements  for  fixed  total  angular  momentum  J  are 


<j  '  £' JMn'|  V  (j£JMn>  =  £  <n'|  VX(R)  jn>fx(  j  '  £' ;  j£;  J)  (5-14 

X 


where  the  fx  are  the  Percival-Seaton  coefficients^® 


£.+£.+J 

fx(j'ft';j*;J)  -  (1)  (jOjO  bjXO)  (£f0£c|  £'£X0)  x  (5-15 

j  '  i '  j 

[(2j'+l)  (2£'+l)  (2 j+1)  (2£+l)]*  {£  j  x  } 

Several  different  potentials  have  been  proposed  to  describe  the  HCl-Ar 
system.  Typically,  they  are  of  the  form  given  in  Equation  5-13.  In  both  our 
calculations  and  those  of  Boulet  and  Robert, the  semiempirical  "A2" 
potential  of  Kircz  et  al.^l  was  used.  In  this  case  X  =  0,  1  and  2  and  the 
radial  potentials  VX(R)  are  6-12  Lennard-Jones  potentials  whose  parameters 
were  obtained  by  fitting  half-widths  calculated  using  semiclassical  S-matrix 


theory  to  experimental  measurements.  These  radial  functions  are  displayed  in 
Figure  5.1. 


5.2.2  Dipole  Autocorrelation  Function 

Since  the  Hamiltonian  is  diagonal  in  the  total  angular  momentum 
(Eq.  5-12),  the  single-perturber  dipole  autocorrelation  function 

$(t)  =  Tr  [y  {elHT/R  y  e_lHT/R  p(H) } ]  (5-18) 

can  be  expressed  as  a  sum  of  terms  corresponding  to  fixed  initial  and  final 
total  angular  momentum.  Specifically,  for  the  jj  ->  jf  transition  in  the 
uncoupled  line  approximation  and  using  the  zero-order  basis  set  defined  in 
5.1.1, 


*  (t)  =  (1/3)  (-1) 
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For  comparison  with  the  results  of  Boulet  and  Robert, it  is  necessary 
to  relate  the  single  perturber  correlation  function  $  (t)  to  the 
autocorrelation  function  for  a  dilute  mixture  $np(t)  where  np  is  the  perturber 


density.  By  analogy  with  DTC  theory  (Appendix  A),  this  is  accomplished  by 
assuming  that  the  autocorrelation  function  for  Np  perturbers,  $Np(x),  is 
simply  the  product  of  Np  uncorrelated  single-perturber  distributions.  Then, 


N 


Nt 


P  KE.-E  )  x/R  -i(Ef-E.)x/fi 

$N  (x)  =  II  e  $(x)  =  [e  $(x)] 


(5-21) 


where  the  complex  exponential  factor  cancels  out  the  zero-order  time 
dependence  from  Hra( j.  Following  normalization  to  unity  at  t  =  0,  Equation 
5-21  should  be  equal  to  the  normalized  autocorrelation  function  in 
Reference  18,  i.e.. 


or 


-i(E,-E  )x/fi  N_  -n  a  >(x) 
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where  0  is  the  Lennard-Jones  distance  parameter  and  fi  =  Np/np  is  the  system 
volume 


4nR 


Q  = 


(5-24) 


for  a  sphere  of  radius  R0.  In  Equations  5-22  and  5-23,  the  "potential 
correlation  function"  i<i(x)  is  defined  as  the  average  over  velocity  and  impact 
parameter  of  the  time  evolution  of  initial  and  final  zero-order  states. 


As  a  preliminary  test  of  the  RRGM  method,  was  initially  calculated 

by  assuming  the  diagonal  transition  (lj  =  If  and  nf  =  nf  in  Equation  5-19)  are 
the  dominant  contributors  to  the  autocorrelation  function  at  short  times.  The 
results  from  this  calculation  are  shown  in  Figures  5.2  and  5.3  for  times  up  to 
1.0  and  0.4  ps,  respectively.  There  is  a  clear  discrepancy  between  the  RRGM 
and  French  results,  even  for  time  less  than  0.4  ps.  Since  our  analysis  of  the 
convergence  of  the  RRGM  results  with  respect  to  the  number  of  recursions 
indicates  that  our  results  should  be  accurate  (at  least  up  to  about  0.4  ps), 
the  discrepancies  in  Figure  5.3  probably  result  from  a  lack  of  convergence 
with  respect  to  the  size  of  the  basis  set  and/or  the  neglect  of  off-diagonal 
transitions  in  Eq.  (5-22). 

Despite  the  discrepancies,  there  are  also  some  encouraging  features. 

Most  notable  is  the  indication  of  long  time  (impact)  behavior  starting  at 
about  0.6  ps.  The  slope  of  the  approximately  linear  long  time  tail  for  the 
RRGM  results  gives  a  halfwidth  of  0.02  cm-*  which  is  in  reasonably  good 
agreement  with  experimental  measurements.  In  addition,  there  is  some 
indication  that  for  times  greater  than  about  0.2  -  0.3  ps  only  a  relatively 
few  transition  amplitudes  are  contributing  to  the  autocorrelation  function. 
Significantly  more  transitions  appear  to  contribute  to  the  short  time  dynamics 
However,  these  transitions  can  be  treated  quite  accurately  with  just  a  few 
recursions.  As  a  result,  although  a  great  amount  of  CPU  time  was  expended  in 
obtaining  the  RRGM  results,  intelligent  optimization  of  the  FORTRAN  codes  in 
terms  of  which  transitions  and  corresponding  numbers  of  recursions  are  used  to 
compute  the  autocorrelation  function  could  significantly  decrease  the  required 
CPU  time. 


6.0  SUMMARY 


Several  critical  military,  geophysical,  and  energy  related  applications 
require  accurate  predictions  of  long  path  transmission  in  atmospheric  window 
bands.  Accurate  predictions  of  transmission  in  these  regions  requires,  in 
turn,  more  accurate  characterization  of  far  off  resonance  absorption  line 
shapes  than  currently  available.  Accurate  experimental  absorption  data  is 
presently  difficult  to  determine  because  of  large  experimental  uncertainties. 
At  the  same  time,  theoretical  approaches  have  typically  been  limited  by  severe 
computational  (cpu  time  and  storage)  requirements. 

This  report  has  presented  an  analysis  of  two  distinct  approaches  to  the 
problem  of  calculating  far  wing  absorption  for  collision  induced  broadening. 
The  first  approach  uses  standard  time  dependent  perturbation  theory  to 
evaluate  the  time  dependent,  dipole  autocorrelation  function,  followed  by  a 
frequency  Fourier  transform  to  give  the  absorption  cross  section.  Two  methods 
of  this  kind  were  evaluated  which  differed  in  the  limit  that  perturbation 
theory  broke  down.  The  second  approach  evaluates  the  autocorrelation  function 
nonperturbatively  using  the  recursive  residue  generation  method  (RRGM) . 

The  perturbation  approaches  closely  resemble  Anderson's  impact  theory  for 
half-widths  in  that  they  are  strictly  valid  only  for  weak  interactions.  These 
approaches  do  offer  some  improvement  in  principle,  in  that  they  treat  all 
degrees  of  freedom  quantum  mechanically  and  satisfy  the  fluctuation 
dissipation  theorem.  In  addition,  although  the  Davies,  Tipping,  and  Clough 
(DTC)  method  relies  on  ad-hoc,  Anderson-like  cutoff  functions  in  the  frequency 
domain  to  ensure  physically  well  balanced  absorption  cross  sections  in  strong 
collision  regimes,  the  generalized  Leavitt-Korf f  (GLK)  method  is  cut-off  free 
in  the  sense  that  the  autocorrelation  function  remains  well  behaved  as  the 
impact  parameter  approaches  zero. 
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The  perturbation  methods  have  been  used  to  treat  H£0  self  broadening. 

H2O  absorption  is  particularly  important  in  the  10  pm  region  where  absorption 
is  thought  to  arise  from  line  wings  of  the  6.3  pm  and  pure  rotation  bands. 
Experimental  measurements  in  this  region  indicate  that  the  x-f actors  (ratio  of 
observed  absorption  to  that  predicted  from  Lorentzian  line  shapes)  range  from 
0.1  to  0.6  and  show  a  decrease  with  increasing  temperature  between  250  and 
400  K. 

Application  of  the  DTC  method  to  one  of  the  stronger  rotational  lines  of 
H2O,  (4,1,4)  -+  (505),  gave  X-factors  at  1000  cm”1  which  were  approximately  an 
order  of  smaller  than  experimental  x-factors  at  the  same  frequency.  In 
addition,  the  calculated  x-factors  showed  a  strong  positive  temperature 
dependence  for  250  <  T  <  400.  The  DTC  results  were  based  on  multipole 
expansion  up  to  and  including  dipole-quadrupole  interactions.  Higher  order 
interactions  would  increase  the  x-factors  and  may  even  change  their 
temperature  dependence.  However,  since  multipole  expansions  often  converge 
slowly,  it  is  unclear  to  what  extent  these  interactions  would  improve  line 
wing  behavior.  Initial  results  (up  to  400  cm”1)  from  the  GLK  method  tended  to 
differ  from  DTC  x-factors  further  off  resonance.  This  difference  was  not 
large  enough  to  cause  us  to  suspect  any  significant  change  in  the  temperature 
dependence. 

Analysis  of  our  results  suggest  two  conclusions.  First,  a  multipole 
expansion  of  the  anisotropic  electrostatic  interaction  is  probably 
insufficient  to  accurately  characterize  far  wing  absorption  for  self-broadened 
H2O.  Second,  strong  collisions  contribute  to  the  absorption  up  to  1000  cm”1 
off  line  center.  The  calculated  absorption  strengths  can  differ  significantly 
depending  upon  the  approximations  or  cut-off  functions  used  to  ensure 
physically  reasonable  absorption  cross  sections  in  strong  collision  regions. 


Based  on  these  conclusions,  the  RRGM  was  explored  as  a  possible  approach 
to  line  broadening.  RRGM  is  a  nonperturbative  method  which  can  treat  more 
general  interactions  than  the  typical  multipole  expansion.  As  part  of  an 
initial  feasibility  study,  the  method  was  used  in  a  restricted  context  to 
determine  HCl-Ar  line  wings.  In  general  the  results  were  encouraging  but  more 
work  is  required  for  applications  to  systems  of  atmospheric  interest. 


Vi S 
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APPENDIX  A. 


DERIVATION  OF  DAVIES,  TIPPING,  AND  CLOUGH  LINESHAPE  FUNCTION 

The  time  evolution  in  Eq.  (2-4)  cnn  be  treated  using  standard  time 
dependent  perturbation  theory  by  first  partitioning  the  Hamiltonian  into 
zeroth-order  and  perturbative  terms 


H  +  H0  +  V 

(A-l) 

and  then  defining  the  time  evolution  operator  U(z)  by 

-iHz/R  -iH  z/R  . 

e  =  e  o'  U(z) 

(A-2a) 

* 

+ iHz/R  +,  =  iH  z  /R 

e  U  (z)  e  o 

(A-2b) 

The  operator  U(z)  satisfies  the  integral  equation 

U(z)  =  1  -f  V(z’)  U(z')  dz’ 

J0 

(A-3) 

...  .  iH  z/R  ,,  -iH  z/R 

V(z)  =  e  o  V  e  o' 

(A-4) 

subject  to  the  initial  condition  that  U(0)  =  1.  Eq.  (A-3) 

iteratively  to  give  the  usual  perturbative  series,  i.e. 

can  be  solved 

U(z)  =  l  (-  b*fZ  dzj  ...  fZ n_1  [V(Zl)  V(z2) 

...  V(zn)]  dzn 

n=o  •'0  *'0  ^ 

(A-5) 

In  DTC  theory, 3  H0  represents  the  zeroth-order  Hamiltonian  for  the 
isolated  molecules  plus  the  kinetic  energy  operators  for  their  relative  motion 
and  any  isotropic  intermolecular  interactions,  while  the  perturbation  V 


represents  the  anisotropic  interactions  which  do  not  commute  with  H0.  The 
zeroth-order  Hamiltonian  can  then  be  further  decomposed  as 


D 


ra 


* 

K'y  -• 


H  =  HT  +  H 

O  K  O 


(  A— 6  ) 


where  H^  is  the  molecular  Hamiltonian  for  the  radiator  and  Hq  contains  the 

remaining  zeroth-order  terms.  Here  we  will  neglect  vibration-rotational 

interactions  and  take  H^  to  be  the  Hamiltonian  for  a  rigid  rotor  with 
eigenstates 

h£  j  jxm  >  =  E.t  |  jxm  >  (A-7) 

where  j  is  the  total  angular  momentum,  m  is  the  projection  on  the  quantization 
axis  and  t  represents  any  additional  quantum  numbers  needed  to  completely 
specify  the  rotational  state. 

Performing  the  trace  over  radiator  states  in  Eq.  (2-4)  and  using  the 
Wigner-Eckart  theorem  to  evaluate  the  dipole  matrix  elements,  the 
autocorrelation  function  becomes 


Q° 

*(t)  =  o  I  I  I  l  (2jx  +  1)  P(E.)  <  j.tJ  H  II X  >  (A-8) 

P  1  f  1*  f'  1  1 

x  <  xj  II  p  II j *x *  >  eiWf ' i  Cllf 


C(t)  =  (2j.  +  1)"  l  l  l  l  I  (jf,  ®£>  1,  m  I  jf,  1,  j1,  n»i) 

ii' ,ff 1  m.  m,  m!  m'  m 

l  f  i  f 


(jf,  m'f,  1,  m  |  j^,  1,  j’,  mj)  x  Tr  [P(Hq)  (A-9) 

<  i\  I  u(z)  I  iL  Ti  mi  >  <  jf  Tf  mf  I  U+<T)  I  if  Tf  m'f  >1 


p(E  )  = 


p(Ho)  - 


e^i 

Tr[e‘eHR] 

Qr 

eA 

-tJH 
e  o 

Tr[e~BHo] 

s 

(A- 10a) 


(A-lOb) 


are  the  canonical  density  matrices  for  the  zeroth-order  Hamiltonians  Hr  and 
H0  and  Tr  is  the  quantum  mechanical  trace  associated  with  Hc. 

For  sufficiently  low  perturber  densities,  it  has  been  shown  that  many 
spectral  features  can  be  reproduced,  at  least  qualitatively,  by  invoking  the 
binary  collision  approximation.  In  this  approximation,  the  Np-perturber 
autocorrelation  function  is  taken  to  be  the  simple  uncorrelated  product  of 
single  perturber  autocorrelation  functions.  Equation  (A-9)  then  becomes 


C(t)  =  [Q(t)  ]np 
ii'.ff'  ii ' , f f ' 


where 


I 


Q(t)  =  (2j  +  if 1  £  III  l  (jf,  m  ,  1,  m  |  j  ,  1,  j  ,  m  ) 

•  i  I  /  /  I  f  f  X  X  X  XX 

Ii  ,ff  m  m,  m.  m.  m, 

i  f  l  f 


(jf,  1,  m  |  j'{,  1,  jj,  mj)  x  Tr  [p(H*) 


(A-12) 


<  ji  TI  “I  •  VZ)  1  jl  Ti  mi  >  <  jf  Tf  mf  I  Us+(t)  1  jf  Tf  “f  >] 

is  the  autocorrelation  function  for  the  interaction  between  the  radiator  and  a 
single  perturber  (designated  by  the  subscript  or  superscript  s). 

The  single  perturber  autocorrelation  function  given  in  Equation  (4-12) 
includes  line  coupling  effects.  When  the  zeroth-order  absorption  lines  are 
well  separated,  however,  the  calculation  can  be  simplified  considerably  by 
invoking  the  uncoupled  line  approximation.  In  the  present  context,  the 
uncoupled  line  approximation  is 


Q(t)  =  6  6  Q(t) 

ii',ff'  ii'  ff*  if 


(A- 13) 


where  Qff(t)  is  given  by 


1=  (2J  f  1)  l  l  l  l  (j  {  m  f  !,  m  |  j  {  1,  j  l  m  J 


m  m^  m^  m^ 


x  (jf,  m’,  1,  m  |  jf,  1,  j4,  m^)  x  Tr  [p(H®) 


(A- 14) 


<  ji  Ti  “1  I  Us(z)  I  ji  Ti  mi  >  <  jf  Tf  mf  I  (t)  I  jf  Tf  mf  >] 


Either  Qii}ff(t)  or  Qif(t)  is  evaluated  in  DTC  theory  by  substituting  the 
perturbation  series.  Equation  (A-5),  for  U(z)  and  U(t)  and  keeping  only  terms 
up  to  second  order  in  the  interaction.  This  procedure  breaks  down,  however, 
for  strong  collisions  leading  to  values  for  the  autocorrelation  function  which 
are  not  physically  reasonable.  Consequently,  in  applications  it  is  necessary 
to  introduce  cutoff  functions  which  ensure  that  the  correlation  functions  are 
always  well  behaved.  In  general,  there  is  no  formal  basis  for  these 
functions . 


APPENDIX  B 


COUPLING  MATRIX  ELEMENTS  OF  THE  LONG  RANGE  ANISOTROPIC  INTERACTION 

In  this  appendix  we  determine  the  coupling  matrix  elements  of  the  long 
range  anisotropic  interaction  V  between  two  molecules  in  an  uncoupled  basis 
set 


|  a>=  |  a j>  |  a2>  |  £m^>  |  £e> 


(B-l) 


consisting  of  molecular  eigenfunctions  for  the  internal  degrees -of- freedom  for 
each  molecule 


al>  =  I  VlV 


(B-2a) 


1 V  ■  I  »2j2-2> 


(B-2b) 


a  spherical  harmonic  |  £m^>  for  the  relative  angular  momentum,  and  radial 
eigenfunctions  |  £e>  for  the  relative  translational  motion.  Following  a  brief 
discussion  of  the  long  range  electrostatic  interaction  desired  by  Leavitt, 
the  Q  functions  (the  squares  of  matrix  elements  of  V  summed  over  all  magnetic 
quantum  numbers)  which  arise  in  a  perturbative  treatment  of  line  broadening 
are  derived  for  arbitrary  orders  in  a  multipolar  expansion  of  V. 


We  start  with  the  electrostatic  interaction  given  by  Leavitt^  for  two 
charge  distributions  separated  by  a  distance  R 


v  =  l  l  (-1)  ( 


2  ,  i  +  2e2j/0n,+£,+i 


')  /R  1  2 


(B-3) 


£1  *2 


x  [  I  ^l’yl’^2,y2  ^  £  +£  y  +y 


(1)  (2) 

where  Qm^i  and  Q£2,y2  are  multipolar  moment  operators  for  the 
molecular  systems  1  and  2,  respectively,  and  C+  is  the  adjoint  of  the 
irreducible  spherical  tensor  given  by 


ct  „<*>  ‘  i  — Y*  „<*> 
,U  (2(1+1)  *’’y 


(B-4) 


where  Y^yCR)  is  an  ordinary  spherical  harmonic.  The  factor  of  the  form 
(§)  is  the  binomial  coefficient  a!/b!(a-b)!  and  (a,  a,  b,  B  |  a,  b,  d,  6) 
is  the  Clebsch-Gordan  coefficient  for  the  vector  coupling  of  a  and  b  to  give 
the  resultant  d.  If  r^j  and  r2j  are  position  vectors  for  charges  qj  and  qj  in 
charge  distributions  1  and  2,  then  the  moment  operators  are  given  by 


QVyi~  ^  ’i*11  Ce1.v1(ri±) 


(B-5a) 


QV*2  =  ]  CV2(r2^ 


(B-5b) 


As  used  in  Equation  (A-3),  the  multipolar  moment  operators  are  each 
defined  relative  to  the  laboratory  reference  frame.  The  spherical  harmonic 
Yfc,y  then  specifies  the  orientation  of  the  5  vector  relative  to  the 
laboratory  z-axis. 


m 


The  electrostatic  interaction  given  in  Equation  (A-3),  is  invariant  to 
the  simultaneous  rotation  of  all  particles  in  the  system.  Most  earlier 
expressions  for  the  electrostatic  multipolar  interaction^  require  the 
application  of  rotation  matrices  to  arrive  at  an  equivalent  form.  This 
invariance  is  easily  seen  by  noting  that  upon  letting  y  =  yj  +  V2,  the 
summation  over  y^  and  y£  becomes 


1  1  [•••l-I  cl<R>+e  l  t1,e2>n1+ii2,ii)  Qt(*J 

"l  “2  u  1  2’  “l  1,1  2’  1 


The  summation  over  y^  corresponds  to  the  contraction  of  the  multipole  tensors 
Q^l^l  and  Q£2jy2  or<^er  ©nd  ^2  *n  th©  molecular  spaces  1  and  2, 
into  a  tensor  Qdj+fc2.y  of  order  dj  +  ^2  in  both  spaces.  Then  using 


ct  (R)  ■  (-l)y  C.  (R) 
d,y  '  ’  d,-yv  ' 


CB-7) 


the  summation  over  y  corresponds  to  a  scalar  product  to  give  a  rotational 
invariant  in  the  spaces  1,  2  and  R,  i.e.. 


I  l  w-qO;2-") 


(B-8) 


The  matrix  elements  of  V  in  the  uncoupled  basis  set,  Equation  (B-2),  are 


V  =  <a’|  V  (ot> 


(B-9) 


a'a 


t,  u,  +  2e,  ,  ,  -<e  +«  +i) 

“II  (-D  <  2t  )’  <‘V|r 


|de> 


«  i  i  ‘1.»2.*i«2.»i+*2)  <»;  j;  -J  ‘C Kw 

“2  '■ 1 


*  <li  Ji  -I  <42’„2  *2*2*2*  <l'*l  K” 

If  the  Wigner-Eckart  theorem  is  applied  to  the  matrix  elements  of 
Q^l,Pl»  Qa2,^2  “d  cal+^2»l1l+JJ2» 

<vi  jl  mi  I  Qn1,p1  ^ihmiy  * 

*<*;  i;  ll  o'"! Iw 


(B-lOa) 


<Jf2  j2  \,V2  ^2  j2  m2>  =  ^2>m2*£2*y2  I  j2’£2’j2*",2) 


(B-lOb) 


<!i 


B-4 


x  (e,m)l,fc1+£2,-(y1+p2)|  £1,£1+£2>£' ,mp 

x  (j1*m1>a1.Vi1l  J fci» J i*mx) 

x  {i2,m2,l2,Vj  j2,£2,j’,o*) 

where  the  factor  ^£^,£2  *s  independent  of  magnetic  quantum  numbers. 


(B-i: 


For  a  perturbative  treatment  of  line  broadening,  we  are  particularly 
interested  in  the  rotationally  invariant  quantities 


«-I  l.l  l.  I  l.  I V. f 

m.  m.  m_  m.  m. 

1  12  2  £  £ 

which  upon  substituting  Equation  (A-ll)  becomes 


(B-l- 


0  -  X  l  1  I  V,.  I  z.«l 

tj  t2 1;  £•  12  12 s  :■  r2  '2 


’£,£„  V  ll 


(B-l! 


where  m  refers  to  the  set  of  magnetic  quantum  numbers  (mi,  m£,  m^) 


Substituting  Equation  (A-13)  for  B 2. ^£2 »  the  summation  over  m  and  m  is 

Vy2  1 

Z  Z  Z  Z  (-D  C-i)  Uj.v^.vJ  VWVW 
y  y  yi  y2 


(£l,yI* £2,y2  K»il2,£i+£2,yI’+y2) 

x  1  1  (^,m)l,2l+a2>y1+p2|  i,fc1+«.2,£'1m£') 


m£  m£ 


x  Y.  Z  ^2 ,m2 ’ £2 * y J  j2’£2’^2,m2^  ^2,m2’  £2>y2  ^2* ^2’^2’^2^ 


m2  m2 


x  Z  I  ( j1,m1,£1,y1l  j j,mj)  (B-16 


ral  ml 


Using  the  Clebsch-Gordan  symmetry  relations,  the  summations  over  mj  and  mj 
are  easly  evaluated  to  give 


I  1  [...]  = 


ml  °1 


BiB.  111  1111  1 


(B-17a) 


x  (jj_.Bj.jj  -mj|  " 


(2j,+l) 


Similarly,  the  summations  over  n»2,  m2'  and  m£,  m£*  are 


(2j1+1)  i 

I  I  [•••] - 1 - 4  «£  £.«u  u. 

,  [(2£,+l)(2£/+l)l  2*  2  Vy2 

^2  n>2  1  1  1  1 


(B-17b) 


I  1  .  I.--]  = 


(2£'+l) 


[(2fc1+2£2+l)(2£^+2£j+l)]  \+*2,*i+*2  \+VylS 


m£  m£ 


(B-17c) 


so  that  the  summation  over  m  and  m'  reduces  to 


(2j1+l)(2j2+l)(2£,+l) 

(2£1+1)  (2£2+l)  (2£1+2£2+l)  6£j,£’  6£2,£j 


(B-18) 
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yl  y2 


(2j1+l)(2j2+l)(2£,+l) 

(2£1+1)(2£2+1) 


a '  o' 
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Then  substituting  back  into  Equation  (A-15)  and  performing  the  summations  over 
and  t2*»  we  find 


(2j  +l)(2j  +l)(2fc'+l) 

<j  =  i5:  — J - 2 - 1  a.  .  i2 


*1*2 


C2)l1+1)(2£2+1) 


iVt2 


(B-19) 
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APPENDIX  C 


RADIAL  MATRIX  ELEMENTS  OF  R" 


When  the  anisotropic  interaction  potential  is  expressed  as  an  expansion 
in  multipolar  moment  operators  for  the  radiator  and  perturber,  the  coupling 
matrix  elements  will  include  radial  factors  of  the  form 


TX  -  <y  I  R -(*+!) 

Ve',  It  ~  <Y£V  1  R 


I  V  > 

1  £e 


(C-l) 


where  X  is  the  sum  of  the  orders  of  the  multipole  tensor  operators  for  the 
radiator  and  perturber  and  are  radial  eigenfunctions  with  energy. 


In  general,  the  eigenfunctions  and  integral  in  Equation  (B-l)  may  be 

determined  by  numerical  integration  (method  1  below).  However,  the 

computation  time  may  be  reduced  by  using  suitable  approximations  to  the 

2 

eigenfunctions  in  the  limits  k^  >  £(£+1)/Pr  (hard  core  regime)  and 
2 

<  £  (£+1)/Pr  (centripetal  regime)  where  Pg  is  the  radius  of  the 
interaction  region  of  the  isotropic  potential.  The  three  basic  methods  used 
to  evaluate  Equation  (B-l)  are  summarized  below.  This  discussion  follows 
closely  the  description  given  in  the  DTC  final  report,^  but  has  been 
generalized  where  applicable  to  arbitrary  values  of  X. 


In  the  current  source  code  MATRIX,  as  well  as  the  previous  GTE  code 
OVERLAP,  the  eigenfunctions  1,£E(R)  are  evaluated  at  a  series  of  equally 
spaced  grid  points  by  numerically  solving  Schrodinger ' s  equation  using 
Numerov's  method. 22  The  integral  in  Equation  (B-l)  is  then  evaluated  using  a 


simple  trapezoidal  or  Simpson  numerical  quadrature.  Typically,  the  radial 
grid  extends  from  R  =  0  out  to  R  =  Rmax  suc^  that 


,_(Xmin  +  15 


max 


<  tol 


where  tol  is  a  user  specified  tolerance  and  Xm£n  corresponds  to  the  lowest 
order  (longest  range)  interaction,  e.g.,  dipole-dipole  for  two  asymmetric 
rotors . 


Method  2: _ Appxoximate  Eigenfunction  in  the  Hard  Core  Regime 

In  the  limit  R  -»  -,  the  eigenfunctions  have  the  form 

sin  [kR  -  fcir/2  +  6  (k)] 

Tfr (R)  -  - * -  (C-2) 

R  •+  -  R 

where  6g  (k)  is  an  energy  dependent  phase  shift. 

At  energies  sufficiently  large  relative  to  the  centripetal  barrier  (i.e 
kz  >  2(2+1)/Pr),  the  radial  wavefunctions  approach  their  asymptotic 
sinusoidal  form  relatively  close  to  the  classical  turning  point  on  the 
repulsive  wall.  Defining  Ra  as  the  radial  position  at  which  the  asymptotic 
limit  is  approximately  valid  for  a  given  pair  of  wavefunctions.  Equation 
(B-l)  may  be  written  as 


where  the  asymptotic  contribution  to  the  matrix  element  is  given  by 


sinfk'R+e' )sin(kR+91 


R 


X+l 


(04) 


and  0  =  6jj(k)  -  fcu/2 


The  asymptotic  integration  may  be  easily  evaluated  for  arbitrary  X  as  follows 
If  we  define  generalized  sine  and  cosine  integrals 


(k,8,x)  = 
m 

I, l  (k,0,x)  = 
m 

then 
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(C-6) 


(C-7) 


For  m  =  1  and  arbitrary  k  and  0, 


Ij  =  -cos (8)  Ci(kJ?a)  +  sin(0)  si(kKa)  (C-9) 

Ij  =  —cos ( 0 )  si(kRa)  -  sin(0)  Ci(kRa)  (C-10) 

where  si  and  Ci  are  usual  sine  and  cosine  integrals  defined  in  Abramowitz  and 

Stegun^^  as 


Ci(x)  =  -  y  dt 

x 


(C-ll) 


si(x)  =  - 


(C-12) 


For  m  >  1,  the  integrals  in  Equation  (B-5)  and  (B-6)  may  be  evaluated  by 
integrating  by  parts.  For  example,  letting 


r  cos(kR  +0) 
_  - B - 


sin(kR  +0) 


(C-13) 


the  cosine  integral  in  (B-5)  becomes  with  successive  integrations 
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or,  more  generally. 
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By  inspecting  Equation  (B-15)  it  is  seen  that  the  last  term  in  the  series  for 
a  given  m  is 


(-l)m^Z  k  A  ,  I  ^  if  m  is  even,  and 
m- 1  l 


(-1)  2  k  A  ,  I,C  if  o  is  odd. 
m-l  l 


Therefore,  letting 


sS  (*> 


£  ("‘1)  A2n+1  ^m-(2n+l) 


(B-18) 


SN  W  "  l  A2n  <fi-2n 


we  have  for  integer  q  =  1,2,  ... 
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IC  -  SC  .  +  SS  ,  +  (-l)q  k  A,  .Is.  for  m  *  2q 

m  q-1  q-1  '  2q-l  1 

IC  -  SC  ,  +  SS  +  (-l)q  k  A.  I?  for  ■  -  2q+l 

m  q-1  q  2q  1  ^ 


(B-20) 


(B-21) 


In  the  absence  of  an  isotropic  intermolecular  potential  V°(R),  the  finite 
solutions  to  the  Schrodinger  equation  are  the  spherical  Bessel  functions 


jt(p)  -  (2tt/p)*  J£+1/2  (p)  p  «  kR 


(B-22) 


which  are  proportional  to  p*^  for  p  -♦  0  and  have  the  asymptotic  form 


sin  (p  -  1/2  Hit) 

h  (P)  - 

P 


(B-23) 


By  comparing  the  asymptotic  limits  in  Equations  (B-2)  and  (B-23)  it  is  seen 
that  outside  the  interaction  region  the  functions  T^g  and  differ  only 
by  the  phase  shift  6£(k).  In  addition,  since  the  shape  of  the  wavefunction 
is  determined  mainly  by  the  form  of  the  potential  to  the  right  of  the  smallest 
classical  turning  point,  the  functions  f^g  and  j£  are  not  very  different 


when  <  £(£+l)/P^  so  that  6fc(k)  is  small  when  £  >  kPjj.  Thus,  in 
the  centripetal  regime  (low  e,  high  £), 


Ti«  W  *  h  <P> 


(B-2A) 


and 


V.'.u 


\  (k'k) ' 


/ 


dR 


J»'  ♦  m  (k,R)  J«  +  m  M) 


(B-25) 


This  is  a  Weber-Schafheitlin  type  integral  which  for  X  >  -1  and  (£*  +  £  -X  + 
2)  >  0  is  proportional  to  a  hypergeometric  function,  i.e.. 


,£e 


S=1 

Z  2  k 


X-l 


max 


2  r(C)r(i-B) 


2F1  (A,B,C;Z) 


(B-26) 


C-7 
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where  T  is  the  gamma  function,  2?l  is  the  Gauss  hypergeometric  series  and 


A  =  i  (£*  +  l  -  X  +  2) 


B  ■  i  (-  0 (£*—£. )  -  X  +  1) 


C  =  A  +  B  +  X 


(B-27) 


Z  =  (k/k')28 

with  0  =  sign  (k'-k)  =  ±  1  depending  on  whether  k'-k  is  positive  or  negative 
and  kniax  equals  either  k'  or  k  depending  upon  which  is  larger. 


For  k'  *  k  (Z  =  1), 


(B-28) 


^(A.B.C;!)  = 


r(C-A)r(C-B) 


Since  C  *  A  +  B  +  X  with  X  a  positive  Integer,  the  hypergeometric  function  may 
be  evaluated  using  Equation  15.3.11  in  Abramowitz  and  Stegun^  for  k'  t  k, 
giving 


IX 

Ve'.te 


Gzl 

Z  2  K 


s,  -  s. 


r(A+x)r(B+x)r(B)r(i-B) 


(B-29) 


s  *  rcxjrd-x)  i  nfrTUHB+ni  d-z)1 

„.0  nl  r(l-X+n) 


(B-30) 


and 


lliHj 


'Vvf 

:5? 


s  =  (z-dx  £  [ 


n!  (n+X)! 


(1-Z)“  x 


(B-31) 


{ln(l-Z)  -  ^(n+1)  -  <*>(n+X+l)  +  «KA+X+n)  +  <KB+X+n)}  ] 


where  $  is  the  diganma  function.  In  calculations  the  summation  in  Equation 
(B-31)  is  carried  out  until  the  relative  difference  between  successive  values 
for  S£  is  less  than  a  specified  convergence  tolerance. 


Siace  X  is  a  positive  integer,  the  parameters  A,  B,  and  C  are  integers 
divided  by  two.  The  parameters  A  and  C  are  always  positive  and  vary  from 
1  to  ^max  +  1  and  fr°ro  3/2  to  £max  +  3/2  respectively,  where  £max  is  the 
largest  accessible  value  of  £.  In  addition,  the  selection  rules  on  changes  in 
the  relative  angular  momentum,  -  X  v  £'-£  v  X,  require  that  B  varies  from  i-X 
to 
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